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Abstract 

We show that the concentration profiles in the subdifTusive system 
with a membrane, which separates a homogeneous solution from a pure 
solvent at an initial moment, has a general scaling form in the long time 
limit C ~ t^F{S/tP), where 5 is a distance from the membrane surface. 
There is also derived the relation involving the subdiffusion parameters 
of the medium and the membrane with the scaling parameters A and p, 
which are measured experimentally. The relation allows one to extract 
the subdiffusion parameters of the system from experimental data. 

PACS numbers: 02.50.-r, 89.75.Da 



1 Introduction 

Subdiffusive systems can be study using the scaling method, which assumes that 
the concentration of a transported substance C{x,t) is of the form 

C{x,t)=t^F{z), (1) 

where 

(2) 

Xq is an arbitrary choosen point (here we consider the one-dimensional system). 
The parameters A and p are controlled by subdiffusion parameters of the system. 
Since the parameters A and p can be measured experimentally, it is possible to 
extract the values of subdiffusion parameters of the system from experimental 
data. This method was used already for subdiffusion in homogeneous systems ^ 
12] and in the systems with chemical reactions 3^ . Recently the concentration of 
the form was observed experimentally in the membrane subdiffusive system 
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We study the subdiffusive system of two vessels separated by a thick mem- 
brane. At an initial moment one vessel contains a homogeneous solute whereas 



in the second one there is a pure solvent. Such a system has experimentally 
studied [SJ |7j . On the basis of theoretical model we show that the concen- 
trations are given by the scaling functions ^ and We also derive useful 
relations combining the subdiffusion parameters of the membrane and solvent 
with the parameters A and p. 



2 The method 



Let us consider the one-dimensional system with a symmetrical membrane which 
is perpendicular to the x axis where a; = and x = d are the positions of the 
membrane surfaces. In the regions (— cxd,0) and (d, cxd) there is a subdiffusion 
with the subdiffusion parameter a and subdiffusion coefficient inside the 
membrane a subdiffusion occurs with the parameters (3 and Dp. We assume 
that the transport of particles inside the membrane is more hindered than in 
the vessels, so a > /?. 

The concentration is described by the subdiffusion equation with fractional 
Riemann-Liouvillc time derivative [H] 



dC{x, t) 
dt 



D 



gi'" d'^c{x,t) 



where rj = a outside the membrane and rj = (3 inside it. To solve the subdiffusion 
equation we take the following boundary conditions: 



J{0-,t)^J{0+,t), 

C{0-,t) ^ AC(0+,t), 
J{d-,t) = J{d+,t), 
XC{d-,t) = C{d+,t), 



(3) 

(4) 
(5) 
(6) 



where J denotes the subdiffusive flux 8 , the dimensionless parameter A controls 
the permeability of the membrane. We note that the boundary conditions lj3J)- 
© were already studied in the case of normal diffusion in membrane systems 
The initial condition we is chosen as: 



C(2;,0) 



Co, X <0, 
0, 2:>0. 



The solutions of the subdiffusion equations in terms of the Laplace trans- 
forms are 



f -I- 7 — (f — 7)e 



.(I + 7)2 - (I 7)2e-2\MA5;d_ 



, (7) 



Co7 
As 



(1 + jje'V^^/Dfii^) + (1 - ^)g-y/s^/Di,{2d~x) 

(I + 7)2 - (I - 7)2e-2yPv^d 



(8) 
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(1 + 7)2 - (1 - 7)2e-2V^^d 



(9) 

where 7 = X^/U^jUps'^^ ")/2^ t:he indexes 1, M, and 2 are assigned to the 
regions (— oo,0), (0, c?), and (d, 00), respectively. To obtain the inverse Laplace 
transform L^^ , at first we find the series expansion of the functions Q-lO in 



, and next we use the following formula jlO| 



(10) 



where 



F [- 



sin [tt {k^ + v)] T (1 + k^i + v) 
A?! 



a 



a > 0, /i > 0, and the parameter v is not limited. 

Since the membrane used in experiments was not transparent for the laser 
beam and consequently the concentration profiles inside the membrane were 
not known) , we focus in the following on the regions outside the membrane. To 
simplify the calculations we consider the above functions in the long time limit. 
According to the Tauberian theorem, the long time limit corresponds to the 
small values of the parameter s in the Laplace transform. For the subdiffusive 
system with relatively thin membrane described in the paper ^Jj the 'long time' 
is evaluated to be longer than 100 seconds. Let us note that in considered system 
7 ^ 00 when s ^ 0. The functions and © in the limit of small s are 



Ci(x,.s) =Co 



C2{x,s) 



D 



Applying the formula (|10|l to eqs. (|11|) and H12|) . we obtain 



C2{x,t) = 



1 



D 



/3 



CoD 



/3+q/2,q/2 



/3+q/2,q/2 



X — d 



(11) 
(12) 

(13) 
(14) 



To illustrate the functions 113|) and p4|l . there are shown few plots of the 
concentration profiles for different times in Fig.l. Here a = 0.75, (3 = 0.5, 
Da = 1 • lO^'', Dp = 1 ■ lO^"', Co = 1, A = 2, and d = 1 (all quantities are in 
arbitrary units). 

In the study ^ it was found that the concentration profiles in the region 
{d,oo) is given by the function with the argument jSJ. The parameter p 
depends on the properties of subdiffusive medium outside the membrane and it 
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Figure 1: The concentrations C for different times (in arbitrary unit). The 
dotted vertical hues represent the membrane surfaces. 
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is controlled by the subdifFusion parameter a whereas the parameter A depends 
on the both of the subdifTusive parameters a and /3. Comparing the function 
(O with eqs. iP) and we obtain 

/3 = A + p, (15) 

and 

p - a/2. (16) 

Since the parameters A and p are measured experimentally we can get the 
subdifFusion parameter inside the membrane from eq. H15|l . 

3 Final remarks 

We observe that the functions Co — C\[x, t) and C2 (x, t) are of the same scaling 
form and they can be expressed by the function 

t^~'^/'^F_i_pj^a/2^a/2 ( ^ J f n ) ' where 5 is the distance from the membrane sur- 
face. Knowing the scaling form we derive the relations H15() and H16|l . 

Presented procedure allows one to calculate the subdiffusion parameter in- 
side the membrane. Our theoretical result fully coincides with the experimental 
one presented in ^ , where it is reported that the concentration profiles obey the 
scaling function ^ for the system containing the membrane made of PEG 2000 
located in the agarose gel solvent, with p = 0.44 ± 0.01 and A = 0.17 ± 0.01. 
The first parameter was determined by means of the time evolution of near 
membrane layers, the second one was extracted from the time evolution of the 
concentration at the membrane surface x = d, which is C2{d,t) ^ t^~°'/^. Sub- 
stituting these values to the relations (|15|) and H16|l . we obtain f3 = 0.61 ± 0.01 
for the membrane under considerations. 
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